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^ ! We propose how to implement precision neutrino mass spectroscopy using radiative neutrino pair 
CN I emission (RNPE) from a macro-coherent decay of a new form of target state, a large number of 
activated atoms interacting with static condensate field. This method makes it possible to measure 
still undetermined parameters of the neutrino mass matrix, two CP violating Majorana phases, the 
unknown mixing angle and the smallest neutrino mass which could be of order a few meV, determining 
at the same time the Majorana or Dirac nature of masses. The twin process of paired superradiance 
" (PSR) is also discussed. 
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Introduction Neutrinos are still mysterious particles: their absolute mass scale (or the smallest 
neutrino mass), the nature of masses (whether they have Majorana or Dirac type masses), and their 
relation to the leptogenesis theory |T], |2] are not clarified experimentally. Experimental efforts 
to unravel these properties are mainly focused on nuclear targets. Nuclear targets, however, are 
problematic at least in one important aspect, the mismatch of energy scale: the released energy of 
nuclear transition is of order several MeV, and this is far separated from the expected neutrino mass 
range of O[0.1]eV. 

We proposed a few years ago the idea of using atomic targets to overcome this difficulty; RNPE 
from a metastable state |e), |e) — )■ \g) + 'y + i^ii^j- This is an elementary process predicted by 
the ordinary electroweak interaction, and its detection opens a path towards the neutrino mass 
spectroscopy [3] , [1] , by precisely measuring the photon energy spectrum, thereby resolving neutrino 
mass eigenstates ,i = 1,2, 3. 

With smaller released energies of atomic transitions, the atomic decay involving neutrino pair 
emission has a demerit of tiny weak rates, unless a new idea of rate enhancement is taken into 
account. Our enhancement mechanism uses a coherent cooperative effect of a large number of atoms 
interacting with a common field [5], A similar idea goes back to the superradiance (SR for short) 
[7j of a single photon emission, where the decay rate from many atoms is in proportion to n'^V, 
the target number density squared times a coherent volume V, unlike the target number nV in the 
spontaneous decay. 

Atoms in a metastable state |e) may have a lifetime for a long time measurement. If these atoms 
further have a developed coherence, macro-coherent two photon emission, called paired superradiance 
(PSR for short), |e) — )■ jf?) +7 + 7, becomes easily detectable [5], its rate oc n^V, with V a macroscopic 
target volume, unlike the case of usual SR limited by \^ oc the photon wavelength squared. PSR has 
a distinct signature: two photons are back to back emitted and have exactly the same energy. 

We propose in this work to use for the target of RNPE a coherent state of atoms interacting with 
static field condensate (we call this as condensate for simplicity). The condensate is a limiting case 
of multiple soliton solutions, as presented below. Both solitons and condensate are proved stable 
against PSR, but unstable for RNPE. 

PSR, emitting a highly correlated pair of two photons, is interesting from points of application 
such as quantum entanglement. Artificial destruction of solitons and condensate, which can be easily 
reahzed by a sudden application of electric pulse (thus abruptly changing the dielectric constant), 
provides the most efficient mechanism of PSR emission yet to be discovered. If we successfully destroy 
solitons for PSR under complete control, solitons may become qubits for quantum computing. 

On the other hand, creation and subsequent long time control of the condensate removes the 
most serious PSR background for RNPE. We compute macro-coherent RNPE rate oc ti^V of conden- 
sate decay and study sensitivity of spectral rates (spectral shape and event rate) to parameters of 
the neutrino mass matrix, most importantly the fundamental parameter of CP violating Majorana 
phases; the parameter of central importance in explaining the matter- antimatter imbalance of the 
universe. RNPE spectrum shape from the condensate decay is time independent after condensate 
formation and the most unambiguous tool for this process. 

Our method uses laser to trigger RNPE at non-resonant frequencies, which should be a great 
merit since the trigger is not destructive to target atoms. 

The natural unit ^ = c = 1 is used in formulas of this paper. 

Effective atomic Hamiltonian and Maxwell-Bloch equation We consider atoms that 
consist of three levels of energies < eg < e^. The state |e), for example ^D2-state of Ba low lying 
levels, is forbidden to decay to \g) by El transition, while El transitions from \p) to \e) and \g) may 
both be allowed. The important part of Hamiltonian is derived |H], |6] by eliminating time memory 
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effects of \p), ^Pi in Ba, and by making a slowly varying envelope approximation of one field mode 
propagating in a direction. The resulting effective Hamiltonian is restricted to two levels, |e) and 
\g), interacting with field E of frequency u and a definite polarization. The 2x2 matrix elements 
/lab ,o,,b = e,g, are Stark energies; a product of two dipole (El or Ml) transition elements to \p) 
times the electric field squared. Dipole transition elements are related to measurable decay rates 
7pa , a = e, from \p) to \a), thus 

^'-- = T^(72 — TT^' (« = e,^/), (1) 

_ _ 37r(epe + epg) I 'Jpelpg /„s 

We ignored the spin multiplicity factor 2Ja + 1 in the relation rf^^ to '-fab- The final PSR rate formula 
should be multiplied by (2Jp + l)/(2Je + 1) if one includes this multiplicity. 

The equation for the polarization vector R (3 bilinears of amplitudes times the target number 
density n), called the Bloch equation, is derived from the Schrodinger equation, and may be written 
as dtR = \E\'^A4R , where elements of 3 x 3 anti-symmetric matrix A4 are linear combinations of ^ab- 
When this equation is combined with the Maxwell equation, written as {dt + drc)\E\'^ = ijjfige\E\'^R , 
with R a component of R, a closed set of equations follows, to describe spacetime evolution of 
polarization and propagating field |8j, |6]. 

When relaxation processes are ignored, one can introduce the tipping angle 6{x,t) by R{x,t) = 
ncos9{x,t). The Bloch equation is then reduced to a relation of 9{x,t) to the electric field strength; 

|i?(x,t)p = dtO/fi with /i = ^ {flee — fJ'ggY + 4/i^e/4- The field 9{x,t) is an analogue of the area for 

field propagation in the two-level problem [9]. The Maxwell equation in terms of 6{x,t) is 

{dt + d,)e = ami-cose + A), (3) 

V ^pe^pg ^pe + ^pg 

where e^a = ^b—^a is the atomic energy difference. For the Ba D-state, ~ 2.4xl0~^cm~^(?7,/10^^cm~ 
at a; = eeg/2. Both and /i depend on oj. The non-linear equation (|3]) describes dynamics of a 
fictitious pendulum under friction periodically varying oc a^sin^ at its location 9. 

For \A\ < 1, the tipping angle is restricted to a finite ^—region of < 27r. The propagation problem 
in this case has been analytically solved in |6j in terms of arbitrary initial data. Hence the system 
appears integrable in the mathematical sense. Typical solutions describe multiple splitting of pulses 
and their compression when they propagate in a long coherent medium, as fully explained in [B]. 
The number of split pulses is given by the initial pulse area 6'(— oo, oo) divided by 27r. This behavior 
of pulse in medium is a symptom of instability, and pulses stabilize via PSR. It is thus anticipated 
that stable objects against PSR exist; solitons. 

Soliton solutions There are two types of analytic solutions for solitons; \A\ = 1 giving a single 
soliton of quantized area 27r and |y4| > 1 the multiple soliton. The case |y4| < 1 is unphysical since 
an excited state of population ncos6' 7^ —n exists at ,^ = ±00. The case of \A\ = 1 solution of area 
27r has been obtained in [6] by using a different method. 

We look for soliton solutions by assuming one variable dependence of x — vt for a soliton of velocity 
V and by reducing the partial differential equation to an ordinary one. The solution for A = 1 thus 
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obtained has a Lorentzian shape of flux and the population given by 



m.,tt='-^ „ , (5) 

/i a^[x - vt)^ + [I - vy 
cos e{xt) = + (6) 

The sohton size is 0[l/am], and its field fiux is of order, Om/yU ~ 30Wmm~^(n/10^®cm~'^) for the Ba 
soliton at oj = eeg/2. 

This method applied to the |y4| > 1 case, on the other hand, gives a new class of solutions given 

by 

\T?f +M2 V AcOsX-1 

\E{x,t)\ = — -, cos^x,t = — i (7) 

/i 1 — V A — COS X A — COS X 

amVA"^ - 1 

A = [X — ft) . (8) 

Unlike the single peak for \A\ = 1, the field fiux given by ([7]) has infinitely many peaks equally 
spaced, describing multiple soliton solutions in medium. 

For a finite length of medium one may impose the boundary condition of no excited state at two 
target ends of x = ±L/2. This gives a condition, amLy/ — 1/ (1 — f ) = 2n{2ns — 1) , = 1, 2, ■ ■ ■ . 
The quantity \/ A^ — 1 is of order, and the soliton number density ~ ns/{amL). 

Solitons may both emit and absorb photons within medium, their rate difference oc cos 9\E{x^ t) p. 
This quantity, when integrated in the entire medium supporting a soliton, gives an integral of a total 
derivative oc S^jSin^, hence vanishes for the quantized area of A6' = 2ti. This proves the soliton 
stability against PSR. 

Field condensate One may consider the limit of large soliton density, ng/amL ~ yj A^ — l/An — )■ 
oo, simultaneous with the limit f — )■ 0. Denoting A = ijlhY + 1 with 77 kept constant, one has 

l^cMP = TTTT ~ (9) 

^ [t] — V cos[amri[x — vt) / v)) fi 

thus an almost constant field fiux is derived. The population oc cos 6 oscillates, with the time period 
r = 2n/{am'J]) and the space period tv. The parameter 77 is 47rx soliton density x soliton velocity. 

Practically, the shortest spatial period is limited by the inter-atomic distance d. By identifying 
the period tv with d, one finds r] ~ 27iv/{amd), hence t = d/v. As w — )■ 0, r — )■ 00, and the 
target becomes fully excited with cos^^ = 1. For the Ba ^D2-state, the relevant numerical value is 
amd ~ 4 X 10-^(n/10i^cm-3)2/3 _ 

The limit taken here gives a constant field rjam/lJ' and the full excitation of target everywhere 
(strictly, this is true for an infinitely long medium). This is the state of field condensate we use for 
RNPE. Field condensate can be created by trigger laser irradiation from multiple directions, since it 
has no memory of a particular direction. 

The stability analysis around the condensate can be made, taking E = Ec + SE ,6 = 6^ + SO with 
Ec,6c = the condensate solution. By keeping linear terms oc 6E, 66 in the Maxwell-Bloch equation, 
with 6E, 66 oc e~*'^* for time dependence, the perturbation equation dx6E = i{u + a^ri/uj)6E gives 
a bounded and purely oscillatory solution, indicating the stability of field condensate. 

PSR rate at soliton and condensate destruction We first mention PSR rate with- 
out soliton creation. The PSR rate without trigger is /i^geggn^y/(2^7r^), which is numerically ~ 
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0.5MHz(n/10^^cm ^yv/cm^ for Ba. Under a strong trigger of flux l-Ep, the rate for a target of 
length L becomes |B] 



fi%e,yVL\E\'' 



327r 



(10) 



Although the rate for liJp ^ lO^Wcm"^ is large, time structure of PSR is complicated [6]. 

PSR after soliton formation occurs only at its destruction, without absorption from \g). The 
emission rate from |e) is oc (1+cos 6)\Es\'^/2. One may compute rates based on perturbative methods, 
in which one of the photons belongs to the soliton pulse. The other photon is emitted backward to 
the soliton propagation direction, with exactly the same energy. The large rate enhancement oc n'^V 
is understood by the momentum conservation among emitted particles, implying = l. 

The PSR rate at soliton destruction is (taking L = dx in eq. ( ITOl) ) 



IblTjJ, Ctmi^ ~ '^'-J + (l ~" "^J 



The rate remains large during a time of 



1 — V 1 — V lO-'^^cm ^ , , 

At = ~ 14/isec . (12) 

varr, V n 



The space integrated rate per soliton 



IS 



eegVnV vr?V 
16/. ^^lO^W^' 

(numbers for Ba) a formula valid for a target of length L ^ l/ftm- For a short target of L < l/om 
the rate is reduced by amL/ (7r(l —v)). The integrated rate for long target is by many orders (~ 10^) 
larger than the trigger-less PSR rate. The prolonged time of 0[l]/isec and its simple profile structure 
has a number of merits of easier PSR identification such as the back to back two photon coincidence 
measurement. PSR rate at condensate destruction is larger by rjamL/lvTr) than at the single soliton 
destruction. 

Effect of relaxation There are a number of processes that might destroy coherence. One of 
them is given by a field decay, introduced by a term in the Maxwell equation. This modifies the 
basic equation ([3]) by an additional term —k,9. With the ansatz of variable dependence of ^ = x — vt, 
this equation is 

(1 - v)^ = -amCos6 - k9 . (14) 
d^ 

Direct numerical integration of eq. ffT^ gives distorted quasi-soliton solutions. Their profile, al- 
though distorted, is unchanged as they propagate. Quasi-solitons exist only for n < where 
Kc ~ 0.725am; indicating a threshold of dissipation. Calculation gives the PSR rate = (PSR rate at 
pure soliton destruction) xA{9 + sm9)/27i (the difference A to be taken at two target ends). This 
rate is smaller than the one without dissipation, but not very much less, unless k is very close to the 
threshold Kc- The condition of a sizable PSR rate, the relaxation constant k < 0[am], implies that 
K < O[0.07]MRz{n/10^^cm-^) for the Ba target. 

RNPE The effective Hamiltonian for RNPE, 

?£!i:^Mf^^.E, (15) 
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gives the amplitude for a single atom, where 5'e and d are electronic spin and dipole operators. To 
give large matrix elements for these, we consider deexcitation of |e) of the angular momentum J = 2 
or J = to 1 5^) of J = via |p) of J = 1, realized by rare gas and alkhali earth atoms. Six measurable 
constants, Ci/s, given by Cij = U*^Uej — %/2 with U the unitary matrix relating the neutrino flavor 
to the mass eigenstate, contain mixing angles and Majorana CP phases |3], |4]. 

The field operator Vi for the Majorana neutrino is a superposition of annihilation {hi) and creation 
{h\) operator of the same Majorana particle, while for the Dirac neutrino it is a sum of two distinct 
operators; particle annihilation [hi] and anti-particle creation {d\). Thus, the UiVj pair emission 



amplitude for i ^ j has the form, hlh^ 



c 



c 



i'^CijOlOj for the Majorana case, and Cijbjd^ 
for the Dirac case. Condensate RNPE decay rate of field \Ec\ ~ r]am/ is a sum of 6 ViVj pair 
emission; 



(16) 



for the Majorana case and Bi 



|cjjp for the Dirac case, while By 



For i 7^ j, Bij 

for both cases. Factors, and /i, attributed to condensate parameters, involve intermediate \p). 
The state \p) that gives the largest condensate factor may be different from the intermediate state 
that gives the largest RNPE rate. In the Yb case, |p) = 6s6p'^Pi for the largest condensate factor 
and |p) = 6s6p^Pi for the largest RNPE rate. 



Yb RNPE rate: D vs M with CP phases 
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Figure 1: Yb RNPE photon energy spectrum in (11) ~ (33) region from condensate decay. Dirac 
case and 3 Majorana cases of different (a,/3) are plotted; Dirac in blue, (7r/2,0) Majorana in dotted 
or short dashed red, (0, 7r/2) Majorana in broken black, and (vr/4, — vr/4) Majorana in dashed purple. 
Neutrino masses of (mi, m2, ms) = (50, 10, l)meV, and cosine ang les 1/^2, 73/2, \/(r97 are assumed 
for all. The Majorana pair emission rate of (a, /3) = (0,0) below (3,3) neutrino threshold is by 
~ 10~^ smaller than the Dirac rate for these masses. Assumed target parameters are n = lO^^cm"^, 
V = Icm^, and rj = 10^. 



The function Iij{uj) in the formula f|T6|) is given by an energy integral arising from the two neutrino 



6 



phase space. The integral in a symmetric form is given in terms of neutrino energies Ei ,i = 1,2 



POO POO I j^y-'-J 3K \ 

=uj J dE,dE.MEi + B2 + w - ee<,)»(C«(Bi, -^2)) I ^ + I . (17) 



f = 2a^iE,,E^ , Kf = -G'SiE„E^^af(E,,E^ + E,E,-S,„n^n, 



- 8 + -A^^ ^ ' ^^^^ 

_ 1 E\ + El - mj - m] {Ej - El - mj + m]f 

where the boundary region is given by Cij{Ei, E2) > with 

a,{E„ E2) = {El + El -ml- m] - - A{El - ml){El - m]) , (21) 

and 5m = 1 for the Majorana and 5m = for the Dirac case. In this calculation, averaged electron 
spin matrix elements, {{k ■ Se/oj)'^) = 1/12 , (S'f) = 3/4, are used. The resulting spectrum given by 
f ll6p sharply rises at each {ij) threshold, a feature characteristic of 3 particle emission of massless 7 
and nearly massless z/j, Uj, when both the momentum and the energy conservation hold. 
The limiting case of 3 massless neutrinos gives RNPE rate of the condensate decay, 

Ghpg^gri^Vy^m 2tu ^ 9(324 - 540x + 245x2) 

^^ej' ^^""^ 32(l-e,,x/(2e,,))2 " ^ ^ 

The coefficient in front of the function f{2u/eeg) is not rapidly varying with the photon energy u in 
the neutrino threshold regions, and on the average over the photon energy to 0.95eV, ~ 8 x 10~^Hz 
for Yb of n = lO^^cm"^ , V = lcm3 , t] = 10^. 

Experiments for the neutrino spectroscopy are to be performed keeping the macro-coherence of 
the condensate. The initial trigger frequency u < Un for RNPE of Uu = eeg/2 — 2ml/eeg, with mi the 
smallest neutrino mass, is reset each time for measurements of rate and parity violating quantities |1] 
at different 7 energies of the continuous spectrum. The energy resolution of RNPE spectrum is thus 
determined by the precision of trigger frequency u, and not by detected photon energy. This is a key 
element for successful implementation of the precision neutrino mass spectroscopy, which must resolve 
photon energies at the /ieV level or less, since the {ij) threshold rise at Uij = eeg/2 — {mi + mjY/{2eeg) 
is separated only a little from the half energy €£3/2 of dangerous PSR. 

Calculated rates are sensitive to Majorana CP phases a, (3 defined by Ue2 oc , t/es oc e*'^ . 
Rate rises at (12), (13), (23) thresholds are oc (|ci2p sin^ a , |ci3p sin^ (3 , |c23p sin^(Q; — (3) ) (the Dirac 
case given by Icjjp without a,P phase), to be further multiplied by an integrated Majorana inter- 
ference factor of (1 — mimj / E1E2) with Ei neutrino energies. For example, 4 cases of (a,/3) = 
(0, 0), (7r/2, 0), (0, 7r/2), (7r/4, — 7r/4) give 3 large threshold phase factors of (sin^ a, sin^ /3, sin^(a — 
/3) ) = (0, 0, 0), (1, 0, 1), (0, 1, 1), (1/2, 1/2, 1) at (12), (13), (23). Since \ci2? ~ 3 x (Idip, \c22?) un- 
der the given oscillation parameters, the magnitude sin^ a is most important in rate comparison at 
low threshold regions, and a small sin^ a gives a much smaller Majorana rate than the Dirac rate 
below (33) threshold. We know of no other measurable quantity of this high sensitivity to a and /3. 
Within a given range of neutrino parameters, the easiest observable might be the Majorana phase, as 
illustrated in our figures. Our proposed experiment is not sensitive to the other CP phase 5, which 
however may be determined by future neutrino oscillation experiments. Determination of all low 
energy phases, a, /3, 5, is a requisite for a better understanding of leptogenesis [2]. 
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YbRNPErate: 13 angle 
with CP phase 
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Figure 2: Yb RNPE photon energy spectrum for different sin^ 6*13 values, 0.03 in blue, 0.02 in 
dotted red, 0.01 in broken black, and in dashed purple, all for {a, (3) = (7r/2,0). Neutrino masses 
of (ms, 1712, mi) = (50, 10, 5)meV, and cosine angles 1/v^, v^/2 for cos 6*23, cos 6*12 are assumed. 
Assumed target parameters are n = lO^^cm"^, V = Icm'^, and rj = 10^. 
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Yb RNPE rate: ml = l,2,4,6 meV 



0.00025 - 



0.00020 - 



0.00015 



0.00010 



0.00005 




1.0694 



1.0696 



1.0698 



eV 



Figure 3: Yb RNPE photon energy spectrum for different rui values; 1 meV in blue, 2 meV in dotted 
red, 4 meV in broken black, and 6 meV in dashed purple, all for {a, (3) = (7r/2,0). Other neutrino 
masses are constrained by neutrino oscillation experiments, and cosine ang les l/^/2, \/3/2, \/a97 for 
cos 6*23, cos 6'i2, cos 6'i3 are assumed. Assumed target parameters are n = lO^^cm"^, V = Icm^, and 
T] = 10^. 
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Yb RNPE rate: Inverted vs normal hierarchy 
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Figure 4: Yb RNPE rate; case of normal and inverted mass hierarchies for different (a, (3) values; 
normal (0,0) (blue), inverted (0,0) (dotted red), inverted (7r/2,0) (broken black) , inverted (vr/4, — vr/4) 
(dashed purple), assuming mi = 5 meV (other masses constrained by oscillation data) and the same 
mixing as FigQ. n = lO^^cm"^, V = Icm'^, 7] = 10^. 



Distinction of Majorana and Dirac neutrinos is possible by the interference effect of identical Ma- 
jorana fermions [3], giving different rates in the vicinity of thresholds. Rate difference of Majorana 
and Dirac pair emission is larger for larger Majorana CP phases, as illustrated in Fig([T]). Experi- 
mentally, the spectral rate is fitted under an assumption of Majorana or Dirac neutrino and either 
hypothesis is verified by a good quality of fitting. 

One possible serious background against RNPE might be the trigger-less SR due to the achieved 
excellent coherence. This process has a monochromatic spectrum at Ceg/S different from RNPE, 
nevertheless it might become dangerous, destroying the initial state. This can be avoided by choosing 
J = — transition, which forbids single photon emission, complete to any order, hence SR 
altogether. Alkhali earth atoms have level structure of this angular momentum configuration. Yb 
and Hg atoms have levels of a similar nature, giving state candidates of |e) = {6s6p)^Pq ,\g) = 
{Gs'^ySo , \p) = (6s6p)^Pi. Incidentally, two photons emitted by — )■ RSR have perfectly correlated 
polarizations, and may serve as an excellent device of quantum entanglement. 

The calculated Yb — )■ RNPE rate for (a, /3) = (0, 0) averaged over all photon energies is 
~ 3 X 10~^Hz for n = lO^^cm"^ , \^ = Icm^,?] = 10^ (a factor to be better understood) and is 
by ~ 70 larger than the corresponding Xe 2 — rate. When the Yb experiment at each photon 
energy u lasts for a day, its event number becomes O[30] if u is in the energy range of Figdl]). This 
event number is further increased by / if one repeats condensate formation with a cycle time of 
1/ f sec. We show in Fig(IT]) and Fig([2]) the spectral rate for various combinations of CP phases and 
the mixing angle ^13. Sensitivity to neutrino masses, in particular to mi values, is shown in Fig(|3]). 
Determination of mi of a few meV range requires a high statistic data near (11) threshold. Distinction 
of normal (3 neutrino masses given by mi,m2 = mf , ma = v/0.052 + 0.012 mf eV's) vs 

inverted (3 masses of mi,m2 = -y^O.OS^ + mf, = a/O.OS^ + 0.01^ -|- ml eV's) hierarchies is most 
dramatic, as seen in Figdl]), hence its determination is easier. 

RNPE rate of condensate decay increases like oc n^, effective with am oc n, as the density n 
increases. The event number from a single soliton decay is smaller than from the condensate decay, 
typically by 1/ {rjamL) for a target length L ^ 1/ am that contains a soliton. 
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In an ideally coherent medium, field condensate never emits PSR. In practice, there may be a 
variety of environmental effects that cause a leakage PSR, a potential background to RNPE. One of 
these effects is a random fluctuation of dielectric constant, most simply due to a density fluctuation 
\f5ri?. The resulting leakage PSR rate is estimated as 



(23) 



for a target length L. We used a Gaussian frequency distribution of width for the trigger. We 
have computed the Yb leakage PSR rate using v6n?/n =5% and Am = IGHz. The calculated Yb 
RNPE near (12) threshold of mi = 2meV is found much larger than the background PSR. The 
leakage PSR becomes larger than RNPE, only at photon energies < 2/ieV away from the first (11) 
threshold. 

In summary, our proposed method of precision neutrino mass spectroscopy is most sensitive to 
Majorana/Dirac distinction and to a, P measurements. It is worthwhile to experimentally investigate 
both formation and long time control of solitons and condensate, which is of crucial importance for 
controlled detection of PSR and RNPE. Some rudimentary method of efficient soliton formation has 
been suggested in [6j. 
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